Stable and convergent fully discrete interior-exterior coupling of Maxwell's equations

Maxwell’s equations are considered with transparent boundary conditions, with initial conditions and inhomogeneity having support in a bounded, not necessarily convex three-dimensional domain or in a collection of such domains. The proposed computational scheme only involves the interior domain and its boundary. The transparent boundary conditions are imposed via a time-dependent boundary integral operator that is shown to satisfy a coercivity property. The stability of the numerical method relies on this coercivity. The method proposed here uses a discontinuous Galerkin method and the leapfrog scheme in the interior and is coupled to boundary elements and convolution quadrature on the boundary. The method is explicit in the interior and implicit on the boundary. Stability and convergence of the spatial semidiscretization are proven, and with a computationally simple stabilization term, this is also shown for the full discretization.